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A THEORETICAL LIGHT-CURVE MODEL FOR THE 1999 OUTBURST OF U SCORPII 
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ABSTRACT 

A theoretical light curve for the 1999 outburst of U Scorpii is presented in order to obtain various 
physical parameters of the recurrent nova. Our U Sco model consists of a very massive white dwarf 
(WD) with an accretion disk and a lobe-filling, slightly evolved, main-sequence star (MS). The model 
includes a reflection effect by the companion and the accretion disk together with a shadowing effect on 
the companion by the accretion disk. The early visual light curve (with a linear phase of t ~ 1 — 15 
days after maximum) is well reproduced by a thermonuclear runaway model on a very massive WD close 
to the Chandrasekhar limit (Mwd = 1-37 ± 0.01 Mq), in which optically thick winds blowing from the 
WD play a key role in determining the nova duration. The ensuing plateau phase (t ~ 15 — 30 days) 
is also reproduced by the combination of a slightly irradiated MS and a fully irradiated flaring-up disk 
with a radius ~ 1.4 times the Roche lobe size. The cooling phase (t ~ 30 — 40 days) is consistent with 
a low hydrogen content X w 0.05 of the envelope for the 1.37M Q WD. The best fit parameters are 
the WD mass Mwd ~ 1-37 Mq, the companion mass Mms ~ L5M (0.8 — 2.0M Q is acceptable), the 
inclination angle of the orbit i ~ 80°, and the flaring-up edge, the vertical height of which is ~ 0.30 times 
the accretion disk radius. The duration of the strong wind phase (t ~ — 17 days) is very consistent 
with the BeppoSAX supersoft X-ray detection at t ~ 19 — 20 days because supersoft X-rays are self- 
absorbed by the massive wind. The envelope mass at the peak is estimated to be ~ 3 x 1O _6 M0, which 
is indicating an average mass accretion rate ~ 2.5 x 1O -7 M0 yr _1 during the quiescent phase between 
1987 and 1999. These quantities are exactly the same as those predicted in a new progenitor model of 
Type la supernovae. 

Subject headings: accretion, accretion disks — binaries: close — novae, cataclysmic variables — stars: 
individual (U Scorpii) 



1. INTRODUCTION 

U Scorpii is one of the well observed recurrent novae, 
characterized by the shortest recurrence period (~ 8 yr), 
the fastest decline of its light curve (0.6 mag day -1 ), and 
its e xtremely helium-rich e jecta (He/H~ 2 by number; see, 
e.g., Webbink et al. 1987 for summary). Historically, the 
outbursts of U Sco were observed in 1863, 1906, 1936, 1979, 
1987, and the latest in 1999. Especially, the 1999 outburst 
was well observed from the rising phase to the cooling 
phase by ma ny observers (e.g., Munari et al. 199E ) , includ- 



ing eclipses ( Matsumoto fc Kato 1999 ), thus providing us 
with a unique opportunity to construct a comprehensive 
model of U Sco during the outburst. 

Our purpose in this Letter is to construct a detailed the- 
oretical model for U Sco based on our light-curve analysis. 
A part of the method has been described in Hachisu & 
Kato (1999), in which they explain the second peak of T 
CrB outbursts, but we again briefly summarize it in §2, 
and will fully describe it in a separate paper. In §3, by 
directly fitting our theoretical light curve to the observa- 
tional data, we derive various physical parameters of U 
Sco. Discussions follow in §4, especially in relation to the 



recently proposed progenitor model of Type la supernovae 
(SNe la). 

2. THEORETICAL LIGHT CURVES 

Schaefer (1990) and Schaefer & Ringwald (1995) ob- 
served eclipses of U Sco in the quiescent phase and de- 
termined the orbital period (P = 1.23056 days) and the 
Ephemeris (HJD 2,451,235.777 + 1.23056.E) at the epoch 
of mideclipse. It is very likely that the companion is a 
slightly evolved main-sequence star (MS) that expands to 
fill its Roche lobe after a large part of the central hydrogen 
is consumed. The phase duration (A<p ~ 0.1) of the pri- 
mary eclipses in quie scent phase in dicates an inclination 
angle of i ~ 80° (e.g., |Warner 1995"| ). 

Our model is shown graphically in Figure l[ We have 
revised Kato's (1990) U Sco light-curve models in the fol- 
lowing two ways: (1) the opacity has been changed from 
the Los Alamos opacity (Cox, King, fc Tabor 1973) to the 
OPAL opacity ( Iglesias fc Rogers 1996| ) and (2) the reflec- 
tion/irradiation effects both of the companion star and the 
accretion disk are introduced in order to follow the light 
curve during the entire phase of the outburst. The visual 
light curve is calculated from three components of the sys- 
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tem: (1) the white dwarf (WD) photosphere, (2) the MS 
photosphere, and (3) the accretion disk surface. 

2.1. Decay Phase of Novae 

In the thermonuclear runaway model (e.g., Starrfield, 
Sparks, & Shaviv 1988), WD envelopes quickly expand to 
~ 10 — 100 i? Q or more and then the photospheric radius 
gradually shrinks to the original size of the white dwarfs. 
Correspondingly, the optical luminosity reaches its maxi- 



summed up the contribution of each patch, but, for sim- 
plicity, we neglect both the limb-darkening effect and the 
gravity-darkening effect of the companion star. Here we 
assume that 50% of the absorbed energy is reemitted from 
the companion surface with a blackbody spectrum at a 
local temperature. The original (nonirradiated) photo- 
spheric temperature of the companion star is assumed to 
be T ph.M.s = 4600 K becau se of B — V = 1.0 in eclipse 



mum lat the maximum expansion of the photosphere and 
then decays toward the level in the quiescent phase, keep- 
ing the bolometric luminosity almost constant. Since the 
WD envelope reaches a steady-state after maximum, we 
are able to follow the development by a unique sequence of 
steady-state solutions as shown by Kato & Hachisu (1994). 
Optically thick winds, blowing from the WD in the decay 
phase of novae, play a key role in determining the nova du- 
ration because a large part of the envelope mass is quickly 
blown off in the wind. 

In the decay phase, the envelope structure at any given 
time is specified by a unique solution. The envelope mass 
AM is decreasing because of both the wind mass loss at a 
rate of M w ; n d and the hydrogen shell burning at a rate of 
M n „n, i.e., 



dt 



-AM = M a , 



M wi 



wind 



M n 



(1) 



where M acc is the mass accretion rate of the WD. By in- 
tegrating equation (|l|), we follow the development of the 
envelope and then obtain the evolutionary changes of phys- 
ical quantities such as the photospheric temperature T p h, 
the photospheric radius R p h, the photospheric velocity 
u p h, the wind mass loss rate M w ; n d, and the nuclear burn- 
ing rate M nuc . When the envelope mass decreases to below 
the critical mass, the wind stops, and after that, the enve- 
lope mass is decreased only by nuclear burning. When the 
envelope mass decreases further, hydrogen shell-burning 
disappears, and the WD enters a cooling phase. 

2.2. White Dwarf Photosphere 

After the optical peak, the photosphere shrinks with 
the decreasing envelope mass mainly because of the wind 
mass loss. Then the photospheric temperature increases 
and the visual light decreases because the main emitting 
region moves blueward (to UV then to soft X-ray). We 
have assumed a blackbody photosphere of the white dwarf 
envelope and estimated visual magnitude of the WD pho- 
tosphere with a window function given by Allen (1973). 
The photospheric surface is divided into 32 pieces in the 
latitudinal angle and into 64 pieces in the longitudinal an- 
gle as shown in Figure [j]. For simplicity, we do not consider 
the limb-darkening effect. 

2.3. Companion's Irradiated Photosphere 

To construct a light curve, we have also included the 
contribution of the companion star irradiated by the WD 
photosphere. The companion star is assumed to be syn- 
chronously rotating on a circular orbit and its surface is 
also assumed to fill the inner critical Roche lobe, as shown 
in Figure [l]. Dividing the latitudinal angle into 32 pieces 
and the longitudinal angle into 64 pieces, we have also 



{ Schacfer fc Ringwald 199E ; see also Johnston fc Kulkarn: 
1992D . Two other cases of r ph ,MS = 5000 and 4400 K have 



been examined but do not show any essential differences 
in the light curves. 

2.4. Accretion Disk Surface 

We have included the luminosity coming from the ac- 
cretion disk irradiated by the WD photosphere when the 
accretion disk reappears several days after the optical max- 
imum, i.e., when the WD photosphere shrinks to smaller 
than the size of the inner critical Roche lobe. Then we 
assume that the radius of the accretion disk is gradually 
increasing/decreasing to 



aRl, 



(2) 



in a few orbital periods, where a is a numerical factor indi- 
cating the size of the accretion disk and R\ is the effective 
radius of the inner critical Roche lobe for the primary WD 
component. The surface of the accretion disk absorbs pho- 
tons and reemits the absorbed energy with a blackbody 
spectrum at a local temperature. Here, we assume that 
50% of the absorbed energy is emitted from the surface 
while the other is carried into the interior of the accretion 
disk and eventually brought into the WD. The original 
temperature of the disk surface is assumed to be constant 
at T ph . disk — 4000 K including the disk rim. The viscous 
heating is neglected because it is much smaller than that 
of the irradiation effects. 

We also assume that the accretion disk is axisymmctric 
and has a thickness given by 



h — /3i?disk 



V -Rdisk 



(3) 



where h is the height of the surface from the equatorial 
plane, w the distance on the equatorial plane from the 
center of the WD, and (3 is a numerical factor showing the 
degree of thickness. We have adopted a ro-squared law 
simply to mimic the effect of the flaring up of the accre- 
tion disk rim (e.g., Schandl, Meyer-Hofmeister, & Meyer 
1997) and the exponent does not affect the disk luminos- 
ity so much mainly because the central part of the disk 
is not seen. The surface of the accretion disk is divided 
into 32 pieces logarithmically and evenly in the radial di- 
rection and into 64 pieces evenly in the azimuthal angle as 
shown in Figure |l|. The outer rim of the accretion disk is 
also divided into 64 pieces in the azimuthal direction and 
8 pieces in the vertical direction by rectangles. 

We have reproduc ed the light cu rves in the quiescent 
phase (-B-magnitude, [Bchaefer 1990 ) by adopting a set of 
parameters such as a = 0.7, — 0.30, a 50% irradia- 
tion efficiency of the accretion disk, and i = 80° for the 
WD luminosity of ~ 1000L© (A B = 0.8 and d = 15 
kpc, or Ab = 2.8 and d — 6 kpc, see below), which 
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is roughly corresponding to the mass accretion rate of 
M acc = 2.5 x 1O~ 7 M yr- 1 for the 1.37M© WD. 



3. RESULTS 

Figure || shows the observational points and our cal- 
culated V- magnitude light curve (solid line). To main- 
tain an accretion disk, we assume a mass accretion rate 

= 1 X 10" 



7 M yr 1 during the outburst. 



To fit the early linear decay phase (t ~ 1 — 10 days after 
maximum), we have calculated a total of 140 ^-magnitude 
light curves for five cases of the WD mass: M\vd = 1-377, 
1.37, 1.36, 1.35, and 1.3 M , with various hydrogen con- 
tent of X = 0.04, 0.05, 0.06, 0.07, 0.08, 0.10, and 0.15 
of the envelope, where the metallicity Z = 0.02 is fixed, 
each for four companion mass of Mms = 0.8, 1.1, 1.5, and 
2.0Af Q . We choose 1.377M as a limiting mass just before 
the SN la explosion in the W7 model (M Ja = 1.378M Q ) of 
Nomoto, Thielemann, & Yokoi (1984). 

We have found that the early 7 day light curve hardly 
depends on the chemical composition and the companion 
mass but is mainly determined by the white dwarf mass. 
This is because (1) the early-phase light curve is deter- 
mined mainly by the WD photosphere (Fig. 0, dotted 
line) and therefore by the wind mass loss rate and (2) 
the optically thick wind is driven by the strong peak of 
OPAL opacity, which is due not to helium or hydrogen 
lines but to iron lines ( Kato fc Hachisu 1994 ). Therefore, 
the determination of the WD mass is almost independent 
of the hydrogen content, the companion mass, or the disk 
configuration. The 1.37Mq light curve is in much better 
agreement with the observations than are the other WD 
masses. 

The distance to U Sco is estimated to be 5.4 — 8.0 kpc, as 
shown in Figure |^, if we assume no absorption (Ay = 0). 
Here, we obtain 5.4 kpc for the fit to the upper bound 
and 8.0 kpc for the fit to the lower bound of th e observa- 
tional po ints. For an absorption of Ay — 0.6 ( Barlow et 
al. 1981 ), we have the distance of 4.1 — 6.1 kpc, and then 
U Sco is located 1.5 — 2.3 kpc above the Galactic plane 
(b = 22°). These ranges of the distance are reasonable 
compared with the o ld estimates o f the distance to U Sco 
such as 14 kpc (e.g., Warner 1995 ). 

To fit the cooling phase (t ~ 30 — 40 days after maxi- 
mum), we must adopt the hydrogen content of X — 0.05 
among X = 0.04, 0.05, 0.06, 0.07, 0.08, 0.10, and 0.15 for 
M\vd = 1-37 M . This is because the hydrogen content 
X is equivalent to the mass of hydrogen burning in the 
envelope. Therefore, X determines the duration of hy- 
drogen shell burning, i.e., the duration of the midplateau 
phase. For X — 0.05, the optically thick wind stops at 
t = 17.5 days, and the steady hydrogen shell-burning ends 
at t — 18.2 days. This duration of the strong wind phase 
is very consistent with the BeppoSAX supers oft X-ray de- 
tection 19 — 20 days after the optical peak (Kahabka et 



al. 1999) because supersoft X-rays are self-absorbed by 



the wind itself during the strong wind phase. It should be 
noted that we do not use this detection of the supersoft 
X-rays to constrain any physical parameters. Hydrogen 
shell-burning begins to decay from t = 18.2 days but still 
continues to supply a part of the luminosity; the rest comes 
from the thermal energy of the hydrogen envelope and the 
hot ash (helium) below the hydrogen layer. The thermal 



energy amounts to several times 10 ergs, which can sup- 
ply a bolomctric luminosity of 10 38 erg s _1 for ten days or 
so until t ~ 30 days, as seen in Figures || and [|. 

The envelope mass at the peak is estimated to be 
AM - 3 x 10- 6 Af Q for M WD = 1.37M , X = 0.05, and 
Z = 0.02; thus, the average mass accretion rate of the 
WD becomes ~ 2.5 x 1O~ 7 M yr -1 in the quiescent phase 
between 1987 and 1999, if no WD matter is dredged up 
into the envelope. Such high mass accretion rates strongly 
indicate that the mass transfer is driven by a thermally 
unstable mass transfer under the constraint that the com- 
panion star is a slightly evo lved main-sequence star (e.g., 
van den Hcuvel et al. 1992) ). 

The wind carries away about 60% of the envelope mass, 
i.e., ~ 1.8 x 1O _6 M , which is much more massive than 
the observational indication of ~ 1 x 1O _7 M in the 
1979 outburst by Williams et al. (1981). The residual, 
~ 1.2 x 1O~ 6 M , can accumulate in the WD. There- 
fore, the WD can grow in mass at an average rate of 
- 1.0 x 1O- 7 M yr" 1 . 

Our model fully reproduces the observational light curve 
if we choose a = 1.4 and /3 = 0.30 during the wind phase 
and a — 1.2 and (3 = 0.35 during the cooling phase for 
1.37M WD + 1.5 M Q MS with i w 80°. Since we ob- 
tain similar light curves for four companion masses, i.e., 
0.8, 1.1, 1.5, and 2.0 M , we show here only the results 
of Mms = 1.5M . It is almost certain that the luminos- 
ity exceeds the Eddington limit during the first day of the 
outburst, because our model cannot reproduce the super 
Eddington phase. 

4. DISCUSSION 

Eclipses in the 1999 outburst were observed by Mat- 
sumoto & Kato (1999). The depth of the eclipses is 
AV ~ 0.5 — 0.8, and it is much shallowe r than that in 
quiescent phase (AB ~ 1.5, [Bchaefer 1990| ). Their obser- 
vation also indicates almost no sign of the reflection effect 
by the companion star. Thus, we are forced to choose 
a relatively large disk radius that exceeds the Roche lobe 
(a ~ 1.4) and a large flaring- up edge of the disk ((3 ~ 0.30). 
If we adopt a size of the accretion disk that is smaller than 
the Roche lobe, on the other hand, the companion star 
occults completely both the accretion disk and the white 
dwarf surface. As a result, we obtain a deep primary min- 
imum as large as 1.5 — 2.0 or more. 

We have tested a total of 140 cases of the set (a, /?), i.e., 
14 cases of a = 0.7 — 2.0 by 0.1 step each for 10 cases of 
= 0.05 - 0.50 by 0.05 step. The best fit light curve is 
obtained for a = 1.4 and (3 — 0.30. Then a large part of 
the light from the WD photosphere is blocked by the disk 
edge, as shown in Figure [|. The large disk size and flaring- 
up edge may be driven by the Kelvin-Hclmholtz instability 
because of the velocity difference between the wind and the 
disk surface. After the optically thick wind stops, photon 
pressure may drive the surface flow on the accretion disk 
(Fukue & Hachiya 1999), and we may expect an effect on 
the accretion disk similar to what the wind has, but it may 
be much weaker. Thus, we may have a smaller radius of 
a = 1.2 but still have (3 = 0.35. A much more detailed 
analysis of the results will be presented elsewhere. 

It has been suggested that U Sco is a progenitor of SNe 
la (e.g., Starrfield, Sparks, & Truran 1985) because its WD 
mass is very close to the critical mass (Mj a = 1.38M ) for 
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the SN la explosion ( Nomoto et aL 1984 ). Hachisu et 
al. (1999) proposed a new evolutionary path to SNe la, in 
which they clarified the reason why the companion star be- 
comes helium-rich even though it is only slightly evolved 
from the zero-age main sequence. A typical progenitor 
of SNe la in their WD+MS model, Mi,wd = 1.37M Q , 

M 2jM s ~ 1-3M , and M 2 ,MS 2 x 10" 7 Af Q yr" 1 just 

before the SN la explosion, has exactly the same system 
as U Sco. Thus, it is very likely that the WD in U Sco is 
now growing in mass at an average rate of ~ 1 x 10~ 7 Mq 
yr _1 toward the critical mass for the SN la explosion and 
will soon explode as an SN la if the WD core consists of 



carbon and oxygen. 
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Fig. 1. — Configuration of our U Sco 1999 outburst model t ~ 9 days after the maximum, i.e., when the WD photosphere shrinks to 
ijph = I.ORq. The cool component (right figure) is a slightly evolved MS (1.5Mq) filling up its inner critical Roche lobe. Only the north 
and south polar areas of the secondary are heated up by the hot component (1.37 Mq WD, left figure) because a large part of the light from 
the hot component is blocked by the flaring-up edge of the accretion disk. Here, the separation is a = 6.87-Rq, and the effective radii of the 
inner critical Roche lobes are R^ = 2.55-Rq, and R^ = R2 = 2.66-Rq, for the primary WD and the secondary MS, respectively. 
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Fig. 2. — Theoretical light curve (solid line) plotted against time (HJD 2,451,000+), together with the observational points (the filled 
circles represent data taken from Matsumoto & Kato 1999, the open circles from the VS-NET archives), for 1.37M0 WD + 1.5Mq MS with 
an inclination angle i = 80°. The brightnesses of the white dwarf photosphere (dotted line) and of the accretion disk surface (dashed line) 
are also added. The distance to U Sco is estimated to be 5.4 — 8.0 kpc for no absorption (Ay = 0) but 4.1 — 6.1 kpc for Ay = 0.6. U Sco was 
detected by BeppoSAX as a supersoft X-ray source (SSS) 19 — 20 days after maximum (Kahabka et al. 1999). The capital letters "S" and 
"P" mean the secondary minimum and the primary minimum, respectively, in the eclipses observed by Matsumoto and Kato (1999). 
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WD + MS + ACDK d=7.5 kpc X=0.05 
1.37 M WD + 1.5 M MS i=80° 
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Fig. 3. — Same as Fig. 2, but for a part of the light curve that takes place during the transition from the wind phase to the cooling phase 
(t ~ 14 — 22 days). When the optically thick wind stops at t = 17.5 days, the photosphere of the white dwarf envelope drastically shrinks 
from ~ O.li?0 to ~ 0.003-Rq within 1 day. This makes a drop Amy ~ 0.3 near t ~ 18 days (HJD 2,451,253), which is followed by a further 
drop Amy ~ 0.4 caused by a reduction of the accretion disk (ACDK) size. The detection of an SSS by BeppoSAX is very consistent with 
the wind duration because supersoft X-rays emerge only after the massive wind stops. 



